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ABSTRACT 


The  transmission  characteristics  of  200  and  500  statute 
mile  ground  wave  propagation  over  land  are  simulated  by 
simple  networks  and  the  responses  of  these  networks  to 
impulse  and  step  functions  are  computed  and  plotted.  The 
simulators  serve  as  regenerative  equalizers  when  placed 
in  the  return  loop  of  stabilized  feedback  amplifiers . 
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1.  Introduction 


It  is  required  to  discriminate  between  electromagnetic  signals 
caused  by  distant  sferics  and  signals  due  to  man-made  pulses,  by  means 
of  their  transient  behavior. 

However,  the  transients  are  modified  by  the  transient  response  of 
the  long  range  wave  propagation. 

The  present  paper  has  a  double  objective: 

1)  To  describe  analytically  and  graphically  the  distortions  and 
blunting  effect  of  propagation  over  representative  distances 
of  dry  land. 

2)  To  devise  equalizing  circuits  that  minimize  the  transient 
distortion. 

The  computation  is  limited  to  ground  wave  propagation.  This  excludes: 

a.  Line-of- sight  radiation  which  is  nearly  undistorted  and,  in 
its  limited  range,  outweighs  other  modes  of  propagation. 

b.  Sky  waves,  that  is  reflections  from  the  ionosphere.  The 
arrival  of  the  first  sky  wave  is  delayed  by  an  interval  which 
is  slightly  affected  by  the  altitude  of  the  source  but  which 
is  mainly  a  function  of  range  and  of  the  altitude  of  the  re¬ 
flecting  layer1.  This  function  is  derived  in  Section  2.2 

Note:  Although  it  is  excluded  from  the  transient  analysis  of  this  paper, 
the  sky  wave  may  be  utilized.  If  the  altitude  of  the  reflecting 
layer  is  known,  the  delay  of  the  first  sky  wave  pulse  is  a  measure 
of  range  and  thus  aids  in  selecting  the  proper  regenerating 
equalizer. 

2.  Analysis 

2.1  Une-of-sight  Exclusion 

The  distance  to  the  horizon  is: 

Dh  ■/  (Re*h,)a-V  4 /  CVh)5-h*  1-1 

where  Re  ■  effective  radius  of  the  earth, 

h j  ■  antenna  height, 
h  ■  source  altitude. 

For  low  antennas  one  may  approximate 

°h  1-2 

where  the  effective  earth  radius  for  electromagnetic  radiation  is  4/3 
the  actual  radius  or  4600  nautical  miles.  Hence 

A  96  J  h  1-3 


1.  See  Reference  1,  Figure  2. 
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when  D  and  h  are  measured  in  nautical  miles.  Assuming  that  man-made 
pulses  requiring  electromagnetic  analysis  do  not  originate  higher  than  2 
n.  M.,  neglect  of  free  space  radiation  is  justified  for  source  distances 

Dj,  >  130  nM 

2.2  Sky  Wave  Utilization  and  Exclusion 

The  difference  between  the  direct  distance  D  over  curved  earth  and 
the  single  reflection  path  to  a  reflecting  layer  at  altitude  h  and  down 
to  the  receiver  is 


AD  ■  2/  h2  ♦  4(R2  ♦  hR)sin2  D-D 
v  4R 


Assuming  that  both  the  ground  wave  and  the  sky  wave  are  propagated 
very  nearly  with  the  speed  of  light^ 


c  ■  162,000  nM/seCj 

2-2 

the  sky  wave  delay  is 

At  *  —  * 

c 

2-3 

In  the  range  between  130  and  S00  nM,  one  may  approximate 

At  A  £  ♦  £  -  ♦) 

2-4 

and  D  «  R  (x+x3-«0,  with  x  » 

2-5 

In  this  equation,  R  is  the  true  Earth  radius  and  R,  D,  and  h  are 
measured  in  nautical  miles. 

Since  it  is  possible  to  measure  h,  the  height  of  the  reflecting 
layer,  by  conventional  local  tests,  one  may  utilize  the  onset  of  the 
first  sky  wave  as  a  measure  of  range,  provided  that  it  is  clearly  dis¬ 
tinguishable  from  the  ground  wave  signal. 

Typical  values  of  sky  wave  delay,  found  either  from  Figure  2  of  Re¬ 
ference  1  or  from  Equation  2-5,  are  tabulated  below: 


Table  1. 

Range 

Height  of 

Sky  wave  delay 

(n.  M.) 

reflecting  layer 

(micro/  sec) 

(n.  M.) 

200 

50 

163 

400 

50 

95 

500 

30 

36 

For  range  computation,  the  observation  time  should  be  200  to  300y  seconds 
to  include  short  range  signals  and  high  altitude  of  the  ionized  layer. 

Of  these,  only  the  initial  35  may  be  useable  for  source  identification, 
if  the  source  is  distant  and  the  ionization  altitude  is  low. 

The  computations  given  in  the  following  sections  indicate  that  this 
time  may  be  sufficient  for  discrimination  between  sferics  and  man-made 
pulses. 


r 


TM-3494  3. 


2.3  Ground  Wave  Transfer  Characteristic 

The  following  analysis  is  based  on  a  paper  by  J,  R.  Johler  of  the 
National  Bureau  of  Standards.  (Reference  1)  This  paper  derives  a 
highly  complicated  equation  for  the  transfer  constant,  as  a  function 
of  the  conductivity  and  dielectric  constant  of  a  level  ground  surface. 

Numerical  solutions  fcr  typical  dry  soil  are  graphed  in  Figures  7 
and  8  of  Reference  1. 

The  curves  for  200  and  for  500  statute  miles  (174  and  435  n.  M.) 
are  used  as  examples  for  detailed  computation. 


2.3.1  200  Mil*  Transfer  Characteristic 


Numerical  values  of  relative  amplitude  and  phase  for  200  mile  range 
are  given  in  Table  2  below. 

Table  2. 


Frequency 

Amplitude  Relative  level 

Phase, 

Simulation 

(kilocycles/sec) 

nanovolts/m 

db 

Johler 

(denominator) 

Radians 

0.1 

0.4 

0 

2.3 

0 

1 

4 

20 

0.2 

0 

3.2 

12 

30 

0.15 

0.04 

10 

40 

40 

0.2 

0.14 

iOO 

250 

58 

1.5 

1.4 

180 

330 

59 

2.3 

2.4 

500 

80 

46 

4.0 

3.9 

1000 

20 

34 

4.5 

4.3 

! . 3. 1 . 1  200  Mile  Simulating  Network 

The  design  of 

a  simple  simulator 

network  for 

this  range  is  carried 

>ut  step  by  step  t 

;o  illustrate  the  technique. 

Johler' s  curve  has  a  low  frequency  asymptotic  gain  slope  of  6  db/ 
octave,  and  a  high  frequency  slope  of  -12db/octave.  This  requires  an 
energy  density  polynomial  of  the  form 

„  _  kw2 


1  ♦  C2  w5 


♦  C,  wl 
4 


♦  Cc  wfc 

o 


A  close  fit  is  obtained  by  making 

e-  in! - 

1  ♦  0.5  n2  +  0.04  n1* 

with  n  *  :  w0  ■  2tr*105  radians/ second. 

To  make  thi^  physically  realizable,  one  must  make 

e  .  kgn _ 

1  +  ajn  -  bn2  -  cjn3 

This  has  the  power  density 

E  -  kn2 

1  ♦  (a2  -  2b')n2  ♦  ({j2  -  2ac)  n1*  ♦  c^n1* 


3  4 


3-2 


3™  3 


3-4 
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By  equating  the  denominators  of  3-2  and  3-4  tan  by  tana  ona  finds 

a  -  1.412 
b  -  0.7S1 
c  •  0.2 


3-5 

3-6 

3-7 


Two  aquivalant  circuits  producing  this  typa  of  polynomial  ara  shown  as 
Figurss  la  and  lb.  Ths  circuit  per  lb,  baing  tha  simplar  ona,  was  salactad. 

The  circuit  has  five  elements  and  since  aquations  3-5  specify  only  three 
parameters,  it  was  arbitrarily  decided  that 


This  determines 


rj  ■  r2  ■  200  ohm.  3-8 


-  0.0144  micro  farad  3-9 
C3  ■  0.025  micro  farad  3-10 
L2  ■  225  micro  henry  3-11 


The  amplitude  characteristic  of  this  simulator  matches  Johler's  curve  to 
within  the  accuracy  of  curve  reading,  that  is  to  within  less  than  1  db. 

The  phase  characteristic  of  the  denominator  alone  agrees  with  Johler's 
curve  between  1  and  1000  kc.  Hie  numerator  contributes  a  fixed  phase  shift 
of  -0.5*  which  does  not  affect  the  group  velocity  characteristic. 

The  low  frequency  phase  reversal  of  Johler's  curve  could  be  simulated 
by  the  type  of  admittance  in  parallel  with  Cj  shown  in  dotted  lines  in 
Figure  la.  It  was  omitted  from  lb  because  it  does  not  noticeably  affect  the 
transient  behavior  during  the  first  50  microseconds. 


2. 3. 1.2  200  Mile  Equalising  Circuit 

It  is  possible  to  design  a  passive  network  that  inverts  and  thus  equa¬ 
lizes  the  200  mil*  transfer  characteristic  in  a  limited  frequency  rang*. 
However,  it  is  simpler  to  make  us*  of  active  circuits  and  to  us*  the  simu¬ 
lator  itself  as  the  "0"  or  loss  part  of  a  stabilized  feedback  loop. 

Input  and  output  of  the  simulator  may  be  buffered  by  cathode  followers 
or  by  equivalent  transistor  circuits.  The  skirt  gain  and  phase  of  the  am¬ 
plifier  are  shaped  to  insure  stability.  If  the  maximum  gain  correction  is 
chosen  as  25db,  the  overall  characteristic  can  be  made  constant  in  gain  and 
phase  from  about  6  kc.  to  1  Me.  This  would  reduce  transient  distortion  to 
about  1  microsecond. 


2. 3. 1.3  Transient  Response  of  200  Mile  Ground  Wave  Propagation 

The  distorting  effect  of  ground  wave  propagation  on  any  given  source 
wav*  form  is  found  as  follows: 

The  simulator  is  substituted  for  the  actual  propagation.  The  Fourier 
transform  of  the  simulator  is  its  response  to  an  impulse.  The  response  to 
a  step  function  is  the  time  integral  of  the  impulse  response.  It  can  also 
be  directly  derived  as  the  Fourier  transform  of  the  network  characteristic 
divided  by  p  ■  jw. 

If  desired,  on*  can  also  derive  the  transient  response  to  a  ramp  function 
by  a  second  integration  of  the  impulse  response,  or  by  a  second  division  by 
p  before  taking  the  Fourier  transform.  It  is  easy  to  approximate  any  given 
source  function  as  a  sum  of  shifted  impulses,  steps  and  ramps  (linear  slopes). 

The  inverse  process,  that  of  reconstructing  a  source  signal  from  the  dis¬ 
torted  signal  received  after  ground  wav*  propagation  over  known  terrain  and 
known  distance,  has  been  discussed  by  Dr,  J.  W.  F.  Ye2. 

2.  See  Reference  2. 
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It  is  quite  laborious,  therefore  it  is  preferable  to  carry  it  out  elec* 
trically,  in  real  tine  by  means  of  the  equaliser  or  inverted  simulator  as 
described  above. 

The  transient  responses  to  impulses  and  to  unit  steps  are  derived  below, 
in  this  section.  They  are  plotted,  as  curves  A  and  B,  respectively,  in 
Figure  4.  The  Fourier  transformations  are  based  on  Campbell  and  Foster's 
tables3. 

Equation  3*3  is  adapted  to  the  form  used  in  reference  3  by  writing 

E  -  JSE - - -  4-1 

p3  ♦  atp2  ♦  a2p  ♦  a3 

If  p  is  measured  in  radians  per  microsecond  one  has 

a}  -  2.36  4-2 

a2  -  2.82  4-3 

a2  -  1.42  4-4 

This  is  factorized  into 

E  -  4-5 

(P  ♦  bj)  (p2  ♦  b2p  ♦  b3) 

with  b.  -  0.81 

b,  -  1.55 
b3  •  1.53 


2. 3. 1.3.1.  Impulse  Response 
One  may  write 

p  1  b, 

(P*b  j)  (pz+b2p«1>3)  "  (p  ♦$)  S*A2  "  (P*bji  [(p*f ■ 

with  8  ■  0.5b2  ■  0.775 

x  b3  -  0.25b22  »  0.964 

The  first  term  of  5-1  is  of  the  form  listed  in  Reference  3,  Table  1, 
Number  448.1,  Its  transform  is 

^Ut)  *  X  5*n*t  exP  (*8t)  (t  in  microseconds) 

The  second  form  of  5-1  must  again  be  factorized  into 

- -  ..  .  -ki  +  k  p*k 

(P*bj)  [(p+8)  2+X2J  pibj-  Tp^ePTT2”  * 

with  ki  -  0.87 

k2  -  0.644  .  ,, 


5-1 

5-2 

5-3 


5-4 


5-5 

5-6 

5-7 


X 


3.  See  reference  3. 
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The  first  term  on  th«  right  hand  side  of  5-5  corresponds  to  Reference  3, 

Table  1,  Nuaber  438,  with  the  transform 

42(t)  "  “ki  exp  S“8 

The  second  term  corresponds  to  Table  1,  Numbers  448.1  and  449.1,  with  the 
transform  k  k  8 

■  [kj  cosXt  ♦  -2*  i-  sinXt]  exp  (-gt)  ,  5-9 

Summing  Zj,  t2  and  Z3  and  dividing  by  k,  one  obtains 

4,..  ■  -exp(-b,t)  ♦  (cosXt  ♦  l*k|-kj^  *inXt]exp(-8t)  5-10 

(t)  kj  X 

substituting  the  numerical  values, 

■  -exp(-0i81t)  ♦  (cos  0.964t  ♦  1.15  sin  0.964t)exp(-0.775t)  5-11 

*  *  j 

This  is  curve  A  of  Figure u4. 


2. 3. 1.3. 2  Step  Response 


Dividing  4-5  by  p  one  obtains 


6-1 


This  has  the  form  of  the  second  term  of  5-1,  with  negative  sign.  By  using 
the  transforms  listed  in  Section  2. 3. 1.3.1  one  finds 

■  exp(-bjt)  -  [cosXt  ♦  sinXt] exp (-Bt) 

kj  X 

and  numerically 

E^j  ■  exp(-0.81t)  -  (cos  0. 964t  -  0.034sin  0.964t)exp(-0.775t) 
This  is  curve  B  of  Figure  4. 

Note:  One  can  verify  that 


Ect)  ■  k/*w4t 

"  ltV)-3FE(«) 

in  agreement  with  the  fact  that  a  unit  step  is  the  time  integral  of  a  unit 
iqmlse. 


6-2 

6-3 


6-4 

6-5 
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2.3.2  500  Mile  Transfer  Characteristic 

Amplitude  end  pheie  of  the  500  mile  tranefer  characteristic  according 
to  Reference  1  are  given  in  Table  3  below. 

TABLE  3. 


Frequency 

Amplitude 

Gain 

Phase 

kc 

nanovolts/m 

db 

radians 

1 

2 

0 

0.2 

10 

20 

20 

0.4 

32 

40 

26 

0.9 

100 

50 

28 

3 

180 

30 

23.5 

4.5 

500 

1.6 

-2 

6 

1000 

0.1 

•  26 

7.5 

2 . 3 . 2 . 1  200  to  500  Mile  Simulating  Network 

It  Is  possible  to  design  a  single  network  simulating  the  entire  500 
mile  transfer  characteristic.  But,  just  as  the  actual  propagation  pro- 
ceeds  mile  by  mile,  it  is  simpler  to  design  a  supplemental  network  simu¬ 
lating  the  difference,  in  db  and  in  radians,  between  the  200  mile  and  the 
500  mile  transfer  characteristics.  This  reduces  the  number  of  elements 
and  the  loss  in  each  partial  network  and  makes  it  possible  to  equalize  in 
steps,  with  each  partial  simulator  in  a  separate  feedback  loop. 

By  comparing  the  characteristics  of  Tables  2  and  3  one  arrives  at  the 
following  table. 

Table  4:  200  to  500  mile  transfer  characteristic. 


Incremental 

Characteristic 

Simulator 

Characteristic 

Frequency 

Loss 

Phase 

Loss 

Phase 

kc 

db 

radians 

db 

radians 

1 

0 

0 

0 

0 

10 

0 

.2 

0 

0 

32 

2 

.3 

1.8 

.2 

100 

8 

1.4 

8.5 

1.6 

180 

15 

2.0 

14.5 

2.1 

500 

28 

2.7 

29 

2.7 

1000 

40 

3 

40 

2.9 

The  last  two  columns  refef  to  the  simulator  derived  below.  The  200  to 
500  mile  loss  has  a  zero  low-frequency  slope  and  a  12  db  per  octave  high 
frequency  slope.  The  phase  difference  is  zero  at  low  frequencies  and  about 
*  at  high  frequencies •  These  requirements  are  met  by  the  networks  per 
figures  3a  and  3b. 

The  network  per  Pigure  3b  was  chosen  because  it  is  simpler.  The  power 
loss  characteristic  corresponding  to  column  4  of  Table  4,  is 


Pfinl 

p  ■ 

p(out) 


1  ♦  5  n2  ♦  n4 


with  n  defined  as  in  3-2 


7 
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This  is  realised  by  the  complex  amplitude  response 

iUfll  -  (1  ♦  2.19jn)  (1  ♦  0.46jn)  7-2 

E(out) 

The  phase  characteristic 

A4  ■  tan-1  2.19n  ♦  tan'1  0.46n  7-3 


2. 3. 2. 2  200  tO  500  Mile  Equalizing  Circuit 

It  is  recommended  to  obtain  incremental  equalization  by  placing  the 
simulator  per  Figure  3b  as  feedback  loss  into  the  loop  of  a  stabilized 
feedback  circuit  similar  to  that  of  Figure  2. 

The  two  feedback  circuits  connected  in  cascade  will  equalize  the  500 
mile  groundwave  transfer  characteristic.  If  the  gain  is  limited  to  25db, 
the  equalized  band  will  extend  from  6  kc  to  400  kc. 


2. 3. 2. 3  Transient  Response  of  500  Mile  Circuit 

The  transient  response  is  found  by  the  combined  Fourier  transform  of 
the  networks  per  Figures  lb  and  3b ,  connected  in  cascade. 


2. 3. 2. 3.1  Impulse  Response 

The  incremental  response  per  7-2  is  brought  into  a  form  suitable  for 
Reference  3  by  writing: 


E^,J  ‘  (P*8^ 


8-1 


with 


g}  ■  0.29 
g2  ■  1.373 


According  to  5-1  and  5-5,  the  200  mile  characteristic  is 

E  (out)  L.  KlP.+.tL+M 
'  P*b  +  P2+b2p*b3 

Combining  8-1  with  8-4  and  dividing  by  k.  the  overall  response  is 

1  +  K2 

e  ■  -l  ^  +  P  ♦  kT 

P  (P+Wi  (P*8lHP+i*)  (P2+b2p+b3)tp2*(gi*g2)p*giga] 

Factorizing  these  terms  and  utilizing  Fourier  integrals  Numbers  452,  448 
and  449  of  Reference  3  one  obtains 

e^ty  ■  3.4  exp(-0.81t)  •  0.51  exp(-0.29t)  -2.01  exp(-1.375t) 

•(0.88  cos0.964t  ♦  0.875  sin0.964t)  exfp  (-0.775t) 

This  is  curve  A  of  Figure  5. 
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2. 3. 2. 3. 2.  Step  Response 

Since  the  step  function  is  the  time  integral  of  the  impulse,,  one  can 
derive  the  step  response  by  integrating  8-6  with  regard  to  time.  One 
finds 


E(t)  »  -4.2  exp(0.81t)  ♦  1.75  exp(-0.29t)  ♦  1.45  exp  (-1.375t) 
♦  (eos0.964t  -  sin0.964t)  exp  (-0.775t) 

This  is  curve  B  of  Figure  5. 
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3.  Summary  and  Conclusions 

It  has  been  shown  that 

1.  The  distance  of  a  signal  source  can  be  computed  from  the  alti¬ 
tude  of  the  ionized  reflection  layer  and  from  the  interval  be- 
between  ground  wave  and  first  sky  wave.  For  distances  between 
150  and  500  milesi  the  sky  wave  delay  is  between  35  and  150 
microseconds. 

2.  The  ground  wave  transfer  characteristic  over  ground  of  known 
conductivity  and  dielectric  constant  can  be  simulated  quite 
closely  by  one  simple  lumped-constant,  minimum-phase 
network,  or  by  two  such  networks  in  cascade. 

3.  When  the  simulator  networks  are  used  as  loss  circuits  in  the 
return  loops  of  stabilized  feedback  networks,  they  produce 
equalizers  that  flatten  the  amplitude  and  phase  transfer  char¬ 
acteristic  from  6  kc  to  above  400  kc. 

4.  The  transient  response  of  the  unequalized  ground  wave  propa¬ 
gation  over  500  miles  to  both  impulse  and  step  functions  is 
practically  finished  in  less  than  15  microseconds. 

From  the  above  one  may  conclude  that  the  onset  of  the  sky  wave  will 
be  recognizable  for  all  but  the  slowest  of  source  signals.  Hence  it  is 
feasible  to  recognize  the  source  range  and  pick  the  correct  amount  of 
equalization  for  a  good  regeneration  of  the  source  signal. 

The  transient  response  to  step  functions  has  no  zero  crossings  in 
the  initial  50  to  100  microseconds.  The  transient  response  to  im¬ 
pulse  signals  has  one  downward  zero  crossing,  2  to  15  microseconds 
after  the  signal  onset.  It  resembles  the  actual  long  range  records. 

This  makes  it  likely  that  the  source  signals  resemble  impulses 
rather  than  step  functions  and  that  the  first  35  microseconds  may 
suffice  to  bring  out  the  differences  between  single  sferics  and  man¬ 
made  pulses. 

The  transient  response  curves  plotted  in  this  report  may  serve  as 
a  reminder  that  some  features  of  records  taken  at  long  range  are  char¬ 
acteristic  of  the  transmission  path  rather  than  of  the  source. 

W.  J.  Albersheim 
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